Permeability of porous materials determined from the Euler characteristic 
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We study the permeability of quasi two-dimensional porous structures of randomly placed over- 
lapping monodisperse circular and elliptical grains. Measurements in microfluidic devices and lattice 
Boltzmann simulations demonstrate that the permeability is determined by the Euler characteristic 
of the conducting phase. We obtain an expression for the permeability that is independent of the 
percolation threshold and shows agreement with experimental and simulated data over a wide range 
of porosities. Our approach suggests that the permeability explicitly depends on the overlapping 
probability of grains rather than their shape. 
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The understanding of liquid flow through porous ma- 
terials is of huge practical importance in many fields of 
research such as soil and material science or biomedical 
applications (for a detailed review see [1 ). In particular 
when the porous matrix has a non-regular structure, as 
typically encountered under natural conditions, it is an 
open question which structural information determines 
its permeability /c, i.e., the rate at which a fluid can flow 
through the matrix at a given pressure difference AP. 
Among various semi-empirical approaches aiming to re- 
late the matrix properties to its capability for fluid flow, 
the Katz-Thomson model is one of the most established 
[2]. It provides a link between k and the conductivity a 
of the fluid phase within a porous material 

k = cl.^{^), (1) 

where c is a constant which depends on the local pore 
geometry, Ic the critical pore diameter and ctq the fluid's 
bulk conductivity. A particular strength of this model is 
the identification of a characteristic length scale Ic which 
is set by the largest diameter of a spherical particle which 
can freely pass the porous medium. Experimental and 
numerical studies have confirmed the applicability of Eq. 
([T]) and in particular the dependence of k on Ic for dif- 
ferent types of natural rocks [SHS]. In general, it is not 
obvious how a/ao is linked to the structural properties 
of the porous matrix. The most successful approach is 
Archie's law 

-=(P^)\ (2) 

C^O \l-(pcj 

where is the porosity, i.e., the volume fraction of the 
liquid phase, the corresponding value at the percola- 
tion threshold (the highest value for ^ below which there 
is no sample-spanning liquid phase) and /i a critical expo- 
nent. Eq. ([2]) has been extensively studied by numerical 



simulations for discrete and continuum percolation mod- 
els. Most authors observe a universal value /i = 1.3 in 
two-dimensional systems [IHS], however, deviations from 
this value have also been found p!QUT2] . Despite its sim- 
plicity, the application of Archie's law to simulations and 
experiments and in particular to natural porous materi- 
als, is limited because (i) Eq. ([2| is only valid close to the 
percolation threshold and even more important by (ii) 
cannot be obtained from a single sample but is only de- 
fined when the underlying formation (diagenetic) process 
of the porous material is known. These limitations may 
explain the partially controversial results obtained for ja 
in the literature. 

To overcome the above problems, a morphological re- 
construction of porous materials based on Boolean mod- 
els, i.e. randomly placed overlapping grains, has been 
proposed p!3HT5] . For sedimentary rock samples, it has 
been demonstrated that their permeability and mechan- 
ical stiffness are very close to those of Boolean models 
for suitably chosen parameters [13 . These parameters 
can be determined from a specific set of morphological 
measures, so-called Minkowski functionals [16 , which al- 
low to characterize the morphology of random geometries 
and have already been successfully employed for a wide 
range of problems, e.g. dewetting patterns of polymer 
films or the large-scale structure of the universe ||T71 [18] . 
So far, it is not clear whether a direct relation between 
the permeability and the Minkowski functionals exists. 
In addition to reducing the computational effort to ob- 
tain such quantities, a direct link between the permeabil- 
ity and morphological measures would allow for a better 
understanding of fluid flow through porous materials. 
In this Letter we study the permeability of quasi two- 
dimensional (q2D) porous structures of randomly placed 
overlapping monodisperse circles and ellipses. We per- 
formed measurements in microfluidic devices and lattice 
Boltzmann simulations which demonstrate that the per- 
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meability is determined by the Euler characteristic (an 
integral geometric measure related to the connectivity of 
the pore space) of the conducting phase. We introduce 
an expression for the permeability that - in contrast to 
Eq. ([2| - is independent of the percolation threshold and 
only depends on morphological properties of an individ- 
ual sample. The expression is found to be in excellent 
agreement with our data over a wide range of porosities. 
For our investigations we used structures created by 
Boolean models with two different types of grains: 
(i) randomly placed overlapping monodisperse circles 
(ROMC) with constant radius r and (ii) randomly placed 
overlapping monodisperse ellipses (ROME) with con- 
stant aspect ratio a/b = 8 with a and b the long and 
short semiaxis, respectively. Both, the position and (in 
case of the ellipses) the orientation of grains is random 
and uniformly distributed. 

In our experiments, microfluidic channels made of poly- 
dimethylsiloxane (PDMS) with dimensions 8 — 10 pm, 
3.5 mm, 10 — 11mm (height /i, width, length) were cre- 
ated by soft lithography [19] . The center of such straight 
channels contained q2D ROMC and ROME structures 
extending over a length L = 3.5 mm. In our samples we 
set r = 30 pm and a = 84 pm. 

As illustrated in Fig. [l^ we apply a hydrostatic pressure 
in the range of AP = 10 — 50 Pa to the sample by con- 
necting two reservoirs to the in- and outlet and varying 
the amount of injected liquid. In order to determine the 
flow properties through the structures, a diluted aqueous 
suspension of 1.3 pm polystyrene tracer particles is in- 
jected into the channel. From the mean particle velocity 
which is measured by digital video microscopy, we obtain 
the averaged fluid velocity v (for details see [20l[2T]). The 
permeability of the structures is then easily obtained by 
application of Darcy's law 



where r] is the viscosity of the fluid and q = v • (j). 

In addition to experimentally measured values of /c, 
we determined the permeability by numerical simulations 
of the fluid flow through the structures using the mas- 
sively parallel lattice Boltzmann application framework 
waLBerla |22j|23]. The experimental structures were 
modeled by a lattice of size 8 x 4000 x 4000. In the simu- 
lations the pressure was set to AP = 50 Pa, comparable 
to the experiments. 
We studied structures with five different porosities cre- 
ated according to the ROMC and ROME algorithm. 
The porosities were chosen to cover a range between 
(j) ^ 0.85 down to the corresponding percolation thresh- 
olds ^ 0.32 (ROMC) and ^ 0.66 (ROME) which 
are estimates from Monte Carlo simulations [24l [25] . 
Each sample is also characterized by its open porosity 
(j)o corresponding to the volume fraction in which fluid 
flow takes place (i.e. the conducting phase excluding re- 
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FIG. 1. (Color online) a) Sketch of the experimental setup to 
measure the permeability through a porous material, b) Illus- 
tration of the morphological quantities to describe a porous 
sample which consists, in this example, oi N — 9 grains (cir- 
cles) leading to 5 connected solid clusters. The porosity (j) 
corresponds to the area fraction of the entire liquid phase 
(hatched + white inclusion) while the open porosity (j)o is only 
defined by the conducting phase (hatched). The perimeter is 
the length of the dashed line and the Euler characteristic is 
X = — 3 (2 liquid phase components — 5 solid clusters) and 
Xo = —4 (the latter neglecting the liquid phase inside the 
cavity formed in one cluster). Exemplarily we show one re- 
alization for the Boolean process for ROMC o) (j) — 0.406, 
(j)o = 0.375, N = 200, Xo = -18 and for ROME d) = 0.805, 
00 = 0.78, = 80 and Xo = —17. The arrows point to the 
pore throats corresponding to the critical pore diameter Dc. 

gions where liquid is trapped in closed cavities, see Fig. 
^p) and the grain density n = N/L'^, where N is the to- 
tal number of circular or elliptical grains and L the linear 
system size. The different morphological quantities rele- 
vant throughout this paper are summarized in Tab. |lj 
For the q2D structures considered here, it is important 
to realize that the critical pore diameter introduced in 
Eq. ^ can be limited by the sample height h. When the 
two-dimensional critical pore diameter Dc, i.e. the size 
of the largest disk that freely penetrates the structure 
is larger than the latter defines the limiting hydrody- 
namic length scale, i.e. Ic = mm{Dc,h). The values of 
Dc for our samples have been numerically computed from 
the Euclidean distance transform [26l |27] of the struc- 
tures and are also shown in Tab. [ij Two typical struc- 
tures and their corresponding critical pores (red arrows) 
are illustrated in Fig. [l]c,d. 

Fig. [2] shows the experimentally (open symbols) and nu- 
merically determined (closed symbols) permeability (nor- 
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TABLE I. Morphological properties of the studied realizations 
of the ROME and ROMC models. Dc and Ic are given in units 
of lattice sites. 
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(b) ROME 



malized by dc) vs the porosity. As expected, both 
sample types exhibit a decreasing permeability when the 
porosity is lowered towards its corresponding percola- 
tion threshold. Our results are consistent with effec- 
tive medium theories for cr/ao which interpolate between 
Eq. ([2| and the dilute limit H |28] (dashed line). When 
rescaling the porosity according to Eq. ([2|, all ROMC 
and ROME data collapse to a single curve demonstrat- 
ing the validity of Archie's law (inset of Fig. [2]) and that 
the transport properties can be described by a master 
curve even far away from (note, that individual data 
points deviate from this curve due to the finite system 
size). However, the biggest drawback of Eq. ([2| is the 
dependence on the percolation threshold which requires 
full knowledge of the diagenetic process. 

In the following, we propose an alternative expression 
for the permeability which does not require knowledge 
about but only depends on the morphological prop- 
erties obtained from the specific sample. As mentioned 
above, morphological properties can be quantified by the 
use of Minkowski functionals, defined by integral geom- 
etry. With respect to scalar additive functionals, it can 
be proven that only a finite number of functionals pro- 
vide independent shape information [29]. In two dimen- 
sions, the Minkowski functionals are (i) the area, (ii) the 
perimeter between the solid and liquid phase and (iii) the 
Euler characteristic x which is the difference between the 
number of connected components of each phase [14 . We 
also distinguish x Xo^ i-e. the Euler characteristic of 
the conducting, i.e. liquid phase, because k is not affected 
by cavities inside solid clusters, which do not contribute 



FIG. 2. (Color online) Experimentally (open symbols) and 
numerically (closed symbols) measured permeability k of 
ROMC (T, V) and ROME structures (A, A) as a function of 
the porosity 0. All data points are normalised by cl^ . The 
dashed lines corresponds to the interpolation formula given in 
[28 . Inset: k vs. the rescaled porosity, according to Eq. ([2]). 
The increasing deviations between experimental and numeri- 
cal data at small permeabilities are a result of discretization 
errors and numerical artifacts which increase close to (\)c- Due 
to finite size effect some are shifted to negative values. 

to the flow. In Fig. [l]3 we provide a simple geometrical 
interpretation of these morphological quantities. 
Motivated by the observation, valid in various stochas- 
tic processes, that the Euler characteristic vanishes near 
the percolation threshold [30l [31] , we propose the follow- 
ing empirical expression for the permeability of a porous 
structure 

For a q2D geometry as considered here, the local pore 
geometry parameter introduced in Eq. ([T]) is c = 1/12 
which corresponds to the geometrical conflnement be- 
tween two parallel plates [32 . From the geometrical in- 
terpretation of Xo given in Fig. [l]3 it follows that 1 — Xo 
(typically referred to as genus) is equal to the total num- 
ber of solid clusters formed by isolated or overlapping 
grains. Consequently the factor (1 — Xo)/N in Eq. Q is 
the fraction of such clusters per grain within the porous 
material. The limiting cases of Eq. ([2| are correctly re- 
produced. For (j) ^ 1 solid grains hardly overlap this 
leads to (1 — Xo) ^ N . In contrast, when (j) (pc-, a 
single percolating solid cluster will form which leads to 
(1 - Xo)/N ^ l/N ^ 0. 

Fig. |3] shows the permeabilities obtained from experi- 
ments (open symbols) and simulations (closed symbols) 
as a function of (1 — Xo)/N . Independent of the shape of 
the grains forming the porous structure, the data points 
which cover more than two decades, are described by a 
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FIG. 3. (Color online) Experimentally (open symbols) 
and numerically (closed symbols) measured permeability of 
ROMC structures (T, V) and ROME structures (A, A) ver- 
sus (1 — Xo)/N. Error bars are smaller than the symbol size. 
The solid line is a fit of Eq. Q with a — 1.27. Inset: k in 
dependence of the open porosity (j)o. 



single straight line with slope a in accordance with Eq. 
Q. Best agreement is found for a = 1.27 ± 0.09 (solid 
line), close to the critical conductivity exponent /i = 1.3 
(Eq. ©). 

The good agreement between our data and Eq. Q sug- 
gests that independent of the shape of single grains, it is 
their aggregation into solid clusters which dominates the 
permeability in porous materials. In contrast to isolated 
grains, where the liquid can easily flow around, extended 
clusters lead to tortuous streamlines and thus reduce the 
sample's permeability. As a result of the different over- 
lapping probabilities, clusters form in ROMC and ROME 
structures at different N. This rationalizes why the frac- 
tion of clusters per grain enters Eq. Q. The impor- 
tance of extended clusters for the permeability is also 
supported by the microscopic flow fields obtained from 
our numerical simulations. In Fig. [4] we compare the ab- 
solute velocity in the horizontal plane averaged over the 
vertical direction for a ROMC and ROME structure with 
(j)o = 0.4 and 0.85, respectively. At the same open poros- 
ity, elliptical grains form larger less compact clusters than 
spherical grains (insets of Fig. [4^,b and c,d). This leads 
to the formation of stagnant zones with vanishing fluid 
velocity in ROME structures and thus to a decreasing 
permeability. This also explains why k of ROME struc- 
tures is always smaller than that of ROMC structures 
with equal (j)o as shown in the inset of Fig. |3] 
We expect that Eq. Q should also hold in three di- 
mensions and for other types of grains, in particular for 
any grains with convex shapes. For porous structures 
where the grain density is unknown or ill-defined (e.g. 
bio-networks) our expression could be applied by using 
a Boolean reconstruction. This method allows one to 




FIG. 4. (Color online) Contour plots of the simulated veloc- 
ity field for two realizations of each process with equal open 
porosities = 0.4 (a,b) and = 0.85. (c,d). The plots show 
the local average velocity magnitude. A logarithmic colorbar, 
that is identical for a,b and c,d respectively, is used for vi- 
sualization. The unit is ]im/s. As illustrated in the insets 
elliptical grains form larger but less compact interconnected 
obstacles which leads to more tortuous streamlines and stag- 
nant zones where the velocity goes to zero. 



derive an appropriate Boolean model (e.g. polydisperse 
mixture of grains) with an effective grain density h from 
the full set of Minkowski functionals, so that Eq. Q can 
be used [33 . We hope that this work will stimulate fur- 
ther research to explore the range of porous materials 
where Eq. Q can be applied. 

In conclusion we have determined the permeability of 
porous q2D structures of overlapping circles and over- 
lapping ellipses. We find quantitative agreement between 
measured and numerical data and with an empirical ex- 
pression or k based on the Euler characteristic of the 
structure. Our expression predicts k over a wide range 
of porosities and does not require knowledge of the perco- 
lation threshold but is related to the fraction of clusters 
formed by overlapping grains. 

We would like to thank Yujie Li and Jakob Mehl for 
inspiring discussions. 
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